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ABSTRACT 

We have undertaken a Parkes ammonia spectral line study, in the lowest two inver- 
sion transitions, of southern massive star formation regions, including young massive 
candidate protostars, with the aim of characterising the earliest stages of massive star 
formation. 138 sources from the submillimctrc continuum emission studies of Hill et al. 
were found to have robust (1,1) detections, including two sources with two velocity 
components, and 102 in the (2,2) transition. 

We determine the ammonia line properties of the sources: linewidth, flux density, 
kinetic temperature, NH3 column density and opacity, and revisit our SED modelling 
procedure to derive the mass for 52 of the sources. By combining the continuum emis- 
sion information with ammonia observations we substantially constrain the physical 
properties of the high-mass clumps. There is clear complementarity between ammonia 
and continuum observations for derivations of physical parameters. 

The MM-only class, identified in the continuum studies of Hill et al, display 
smaller sizes, mass and velocity dispersion and/or turbulence than star-forming 
clumps, suggesting a quiescent prestellar stage and/or the formation of less massive 
stars. 

Key words: line: profiles - stars: formation - stars: fundamental parameters - stars: 
early- type - ISM: molecules - mascrs. 



1 INTRODUCTION 

Massive stars are dynamical and enigmatic powerhouses 
that shape and drive both their local stellar neighbour- 
hood and the ecology and evolution of their host galaxy. 
Despite this heavy influence, the formation and evolution 
of a massive star is not well understood. Particularly per- 
plexing are the earliest evolutionary stages of massive star 
formation (MSF). The difficulty lies in the unambiguous de- 
tection, identification and characterisation of these stages, 
snapshots of which are difficult to obtain as a result of the 
gene ral rarity of candidates and the rapidity of their evolu- 
tion i|Garav fc Lizanolll999l ). 

Whether or not massive stars are scaled-up analogs 
of low-mass stars is still uncertain. Massive stars exert 
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considerable radiative pressure on the surrounding dust 
and gas, which in principle could halt and reverse spher- 
ical infall of the collapsing protostar. Therefore, a 'sim- 
ple' scaled-up version of low mass star formation is in- 
sufficient for massive stars. There are a number of dif- 
ferent approach es to addres s this di lemma as out li ned i n 
the reviews by lEyans et all (12003): iMenten et all (|2005h : 
IZinnecker fc Yorkd (j2007l ); iBeuther et alj (120071 ). 

As their sophistication increases numerical simulations 
in two and three dimensions show that the radiation pres- 
sure problems associated with spheri cal symmetry can be 
overcome (e.g. iKrumhofz et al.l 120090 . It is still not clear 
however, whether the bulk of the mass that finally ends up 
on the massive star co mes from the monolit hic collapse of 
a single dense core (e.g. iMcKee fc Tanll2003l ) or is accreted 
from surrounding lower density gas which is funneled to the 
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centre of a cluster's gravitatio nal potential where the most 
massive stars are forming (e.g lBonnell et al"1ll998T ). 

The natal molecular cloud, from which high-mass stars 
are formed, is expected to be dense, massive and cold, de- 
tectable only at (sub) millimetre wavelengths. Within the 
molecular cloud, core collapse will be triggered. As the col- 
lapsing protostar gains mass, the gravitational energy will 
serve to heat the core and ionise the surrounding material, 
causing an increase in the luminosity of the core. In terms 
of parameter evolution, the initial protostar will be massive, 
cool and of low luminosity. As the core evolves, it will ac- 
cumulate more mass, and the temperature and luminosity 
of the core will increase. This evolution is seen in low mass 
protostars l|Evans et alj|20()5 ). 

In a search for cold cores tha t would mark the earliest 
stages of massive star formation iHill et all (|2005T l under- 
took a (1.2) millimetre SIMBA0 

continuum emission study 
of sources exhibiting signs of methanol maser and/or ra- 
dio continuum emission - both of which have previously 
proven successful t racers of the earliest stages of massive 
star formation (e.g.lMinier et al.ll2000l ; IWilliams et alj|2004 
iPestalozzi et al.ll2005l ). This SIMBA survey revealed a large 
number of millimetre continuum sources (255) devoid of the 
maser and UC Hn sources targeted. Subsequent submillime- 
tre observations of these sources, which were dubbed 'MM- 
only cores', unveiled submillimetre equivalents for each - of- 
ten revealing the multiplicity of individual MM-only sources, 
and con firmed their asso ciation with cold, deeply embedded 
objects (|Hill et al.ll2006h . 

The aforementioned SIMBA survey detected a total of 
405 millimetre continuum sources, which (based on their star 
formati on associa t ion) could be broken into four classes of 
source. IHill et ail l|2005l ) proposed that each of these classes 
of source could represent a different phase of massive star 
evolution, with the MM-only class a possible example of the 
very earliest stages of massive star formation. A caveat how- 
ever, is whether the MM-only cores are currently undergoing 
or will/can support massive star formation. 

In order to determine the characteristic properties of the 
sources in the SIMBA sample, in particul ar the pre v iously 
unknown and unstudied MM-only cores, IHill et all (|2009l ) 
performed spectral energy distribution (SED) modelling of 
a significant fraction of the sample. SED diagrams are useful 
tools from which physical quantities such as the luminosity, 
mass and temperature can be derived (cf. Iwhitnev et al.l 
120031 ; iRobitaille et~aT] 120061 : IHill et allExih If the obser - 
vational data are well constrained (cf. Minier et al.ll2005T ). 
SED modelling can provide useful estimates of each of these 
parameters, for each star-forming core, which allows estima- 
tion of the evolutionary phase of a young massive star. 

The luminosity, mass and temperature of an astronom- 
ical source are fundamental physical attributes that may be 
used to clarify and characterise the nature of a source, possi- 
bly providing insight into their evolutionary status similarly 
as they do for low mass stars i| Andre et al.ll2000h . Indeed, 
the temperature of a source places important physical con- 
straints on the chemical composition of the core, including 
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which chemical species are present in the core, as well as the 
size and types of grains present in the central star-forming 
cores. 

The method employed for the SED fitting was that of 
Bayesian inference, which enabled a statistically probable 
range of suitable values for the luminosity, mass and temper- 
ature, for each source modelled. As SED modelling is heavily 
reliant upon robust data, it was not possible to usefully con- 
strain each of these parameters for a thorough assessment of 
massive star formation scenarios. In the absence of reliable 
far-infrared data, which would serve to constrain the peak 
of the SED and thus parameters resultant from the fit, ad- 
ditional means of constraining the source parameters, such 
as temperature, are necessary. 

Ammonia is an excellent molecular cloud thermometer 
l|Menten et alJ feoOS) from which the rotational temperature 
and ultimately the kinetic temperature of the source can 
be determined. Ammonia is readily detectab le in quiescent 
dark clouds and regions of low luminosity |Ho fc Townesl 
1983) making it perfectly suited to study the birthplace of 
massive stars. The ratio of the hyperfine components of am- 
monia's signature five-fingered spectrum provides the opti- 
cal depth information of the molecular cloud environment. 
As NH3 is a particularly good probe of high density gas, 
and is more resilient to the effects of depletion than other 
high density tracers (e.g. CS), ammonia observations also 
pro vide information pert aining to the density of the core 
fcf. lLongmore et alll2007l . and references therein). 

We have undertaken an ammo nia spectral line survey of 
sources from the SIMBA sample of lHill et all (|2005) in order 
to obtain an independent determination of their tempera- 
ture an d simultaneousl y examine the robustness of our SED 
fitting (|Hill et al.ll2009h . In addition to facilitating and con- 
straining further SED fitting, these ammonia observations 
provide information pertaining to the molecular abundance 
of ammonia in the cores, as well as the physical parameters 
such as linewidth, rest velocity, column density and virial 
mass. With this information, we seek to address formation 
scenarios for massive stars and identify the star formation 
status/phase of the MM-only core. 



2 OBSERVATIONS 

We have undertaken a multi-inversion transition study of 
ammonia, in the low est two inversion transitions, toward 
the SIMBA sample of lHill et all l|2005h . with the aim of de- 
termining characteristic source properties, such as tempera- 
ture and column density. The sample comprises four classes 
of millimetre continuum sources as alluded to in section [T] 
class M are sources with a methanol maser association, class 
R have radio continuum associations, class MR have both 
a maser and radio continuum source, whilst class MM are 
detected solely from their millimetre continuum emission. 

The ammonia observations were undertaken on the Aus- 
tralia Telescope National Facility (ATNF) operated Parkefl 
radio telescope, using the new K-band receiver which op- 
erates from 16 -26 GHz, during three nights spanning 2008, 
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September 17-19. During this period, 244 sources were ob- 
served simultaneously in both NH3 (1,1) and (2,2). 

The K-band receiver, which was partially commissioned 
in September 2008, was used in conjunction with the digital 
filterbank-3 (DFB) backend which contains dual digitisers 
and Compact Array Broadband Backend (CABB jf| proces- 
sorfl. Although the Parkes telescope is 64 m in diameter, 
only the inner 55 m is used for observations at 23 GHz. 

The NH3 (1,1) and (2,2) observations were undertaken 
in position switch mode, using a bandwidth of 128 MHz cen- 
tred at 23708 MHz to enable simult aneous observations o f 
the two transitions. Maser velocities (jPestalozzi et al.ll2005T ) 
were adopted as the rest velocity for each of the sources 
targeted. For those sources devoid of maser associations 
i.e. class R and MM, the velocity of the nearest coincident 
methanol maser was adopted as the rest velocity. We used 
8192 channels to produce an expected velocity resolution of 
0.2km/s. 

The pointing accuracy of the Parkes telescope is typ- 
ically 10-15 arcsec which is smaller than the beamsize of 
58 arcsec. The Tsys ranged from 85 to 132 K during the pe- 
riod of the observations (12-15hr observing shifts), with a 
median of 96 K. 

Typical integration times were 5 minutes per source. 
The rms noise in the spectra ranges from 0.13 to 
0.30Jy/beam for the NH 3 (1,1) transition and 0.13 to 
0.25 Jy/beam for the (2,2) transition with the exception of 
one source (G 294.989-1.720) which had an rms of 0.86 and 
0.87 Jy/beam for each of the transitions, respectively. The 
full set of spectra, including the sample spectra given in Fig- 
ure [TJ are presented in the Appendix. 



3 DATA PROCESSING 

The data were reduced using the ATNF Spectral Anal- 
ysis Package (ASAPjf) and fit using the Continuum and 
Line Analysis Single-dish Software (CLASS jf] package, using 
standard procedures as described in the respective manuals. 



3.1 Data Reduction 

ASAP is a new software package developed by the ATNF to 
reduce single-dish single-pointing spectral lines observations, 
with particular application to the ATNF suite of telescopes. 

During the reduction procedure, the data were read into 
ASAP, an auto_quotienl[3| was applied, as well as a first or- 
der polynomial before the data were frequency- aligned. The 
data were then averaged over time and polarisation, the gain 



3 The CABB hardware was originally developed for 
the Australia Telescope Compact Array (ATCA) and 
was also used in the DFB3 backend at Parkes. See 
http: / /www. narrabri.atnf.csiro.au/observing/CABB. html 

sec http:// www. parkes .atnf.csiro.au/ observing/documentation / - 
software/ C O RREL/index . html 

5 http:/ /www. atnf.csiro.au/computing/softwarc/asap/ 

6 http://www.iram.fr/IRAMFR/GILDAS/doc/html/class- 
ht ml/class . html 

7 The ASAP auto_quotient command combines the nearest source 
and reference spectra, dividing and subtracting the latter off. 



elevation correction was applied, as was an opacity correc- 
tion of 10 per cent (priv. communication J. Reynolds). The 
rest frequency was then set (23.694512 and 23.7226336 GHz 
for the NH3(1,1) and NH 3 (2,2) transitions, respectively) 
and a third order polynomial was applied. 

Comparison of the main hyperfine components Vlsr of 
the NH3(1,1) and NH 3 (2,2) transitions reveals a small ve- 
locity shift between the (1,1) and (2,2) spectra. The mean 
offset in velocity between the two transitions is 0.24 km/s 
with a standard deviation of 0.26 km/s. As mentioned in 
section[2] the expected resolution of these ammonia observa- 
tions is 0.2 km/s. The rest frequencies, which were obtained 
from the Lovas catalogue, are believed to be less accurate 
than the velocity shift seen in our spectra. 

3.2 NH 3 Fitting 

Each of the NH3(1,1) and NH 3 (2,2) spectra were fit using 
the standard NHa( 1,1) andNH 3 (2,2) methods, respectively, 
in CLASS. Under the assumption that each of the hyperfine 
components have equal excitation temperature and the same 
Gaussian velocity structure, these methods return the calcu- 
lated flux density, rest velocity (Vlsr), linewidth (AV), & 
optical depth (r) as well as the corresponding uncertainties 
to the fits for each source. A 3-<r cutoff was used to define 
non-detections for both transitions. 

While the data quality is mostly of a high standard, a 
subset of the spectra show approximately sinusoidal varia- 
tions in the baseline, typically with period 40— lOOkms - . 
Attempting to use polynomial fitting to remove baseline 
variations with a similar velocity scale as the hyperfine struc- 
ture (~50kms _1 ) may substantially alter the structure of 
the spectra. While the line Vlsr and AV are likely robust 
(the line width of the main component and satellites are 
typically much smaller than the baseline variations) , the ra- 
tios of measured flux density, both for the hyperfines within 
a transition and between different transitions, are not. De- 
rived optical depths and temperatures should therefore be 
treated with caution for these sources. Any spectra showing 
this anomaly are flagged as BR (baseline ripple) in column 
11 of Table [T] Some sources, particularly those towards the 
Galactic centre, show multiple blended velocity components, 
making the fits less reliable. These data are flagged as BL 
(Blended) in column 11 of Table [1] All other spectra are 
considered reliable and marked as an R in this column. 

The parameters derived from the fits to both the 
NH 3 (1,1) and NH 3 (2,2) spectra are given in Tabled] The 
source name is given in colu mn 1, in ri g ht asc ension order, 
using names consistent with iHill et all (|2005l ) . The source 
class is given in column 2. Columns 3, 4, 5 and 6 present the 
flux density, Vlsr, linewidth (AV) and the optical depth 
for the NH 3 (1,1) transition. Columns 7-10 present these 
parameters for the NH 3 (2,2) transition. For sources with an 
NH 3 (1,1) detection but no NH 3 (2,2) detection, a 3-a upper 
limit is derived from the NH 3 (2,2) spectra which is given in 
column 7, and no values are listed in columns 8-10. 

Column densities and rota tional temperatures we re de- 
rived in the standard way (e.g. lUngerechts et al.lll986T ). The 
kinetic temperature estimates were calculated from the ro- 
tational tem p eratu re following the procedure outlined in 
iTafalla et al.l ([2004!) . The derived molecular gas properties 
for each source are presented in Table [2] in right ascen- 
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Table 1. NH 3 (1,1) and (2,2) fit results. 
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Figure 1. Sample NH3(1,1) and (2,2) spectra, shown left to right, respectively. The y-axis is in units of Jy/beam. The green lines on 
the spectra indicate the fits as determined with CLASS. Top Row: G9. 88-0. 75 a class M (see section[2} source, and bottom: G8. 718-0. 410 
a class MM source. The full set of NH3(1,1) and (2,2) spectra can be found in the Appendix. Whilst a sample of the full appendix is 
included here (Fig. ??), the entire appendix of spectra can be found in the online version. 



sion order. The rotational temperature (T rot ) and kinetic 
temperature (Tkin) are given in columns 2 and 5, with 
the lower and upper limits to each temperature given in 
columns 3 & 4 and 6 & 7, respectively. Temperature uncer- 
tainties are discussed further below. The column density of 
the NH3(1,1)(Nj v jjj 3 ^ js) and the total column density of 

the gas (Njyjjj^) are presented in columns 8 and 9, respec- 
tively. 

Given the varying robustness of the NTi3(l,l) and (2,2) 
detections, sources were divided into 4 groups based on 
the reliability of the spectra (and hence derived physical 
properties). Group 4 sources were defined as those with no 
NH3(1,1) detection. No physical properties were derived for 
these sources, which can be found in Table[3] Sources with an 
NH 3 (1,1) detection, but with either i) no NH 3 (2,2) detec- 
tion or ii) the NH 3 (1,1) & (2,2) Vlsr/AV differed by >3 
km/s (making it unlikely the emission is coming from the 
same gas) were defined as group 3. For these sources, there 
is no Nfi3(2,2) emission at the position of the NH3(1,1) de- 
tection, so the RMS of the NH3(2,2) spectra was used to 
derive an upper temperature limit and column density esti- 
mate. These numbers are highly uncertain. Group 2 sources 
are defined as those with both NH3(1,1) and (2,2) detections 
which have consistent kinematics but the signal to noise of 
at least one spectra lies in the range 3 < a < 10. Finally, 
sources with both NH3(1,1) and (2,2) detections >10<7 are 
defined as group 1. These group allocations are given in col- 
umn 10 of Tabled 

Uncertainties in the temperature were estimated for 
group 1 and 2 sources using the uncertainties in the mea- 



sured parameters (Flux density, Vlsr, AV & r) from the fits 
to the spectra. For each core, a maximum and minimum ro- 
tational temperature were derived assuming each of the pa- 
rameters was at la above or below the actual measured value 
(e.g. Flux 11 TAFlux 11 or Flux 11 -AFlux 11 ). The lower and 
upper rotational temperatures are shown in columns 3 and 
4 of Table [2] These extrema were then used to derive cor- 
responding lower and upper limits t o the kinetic tempera- 
tures, again following iTafalla et al.l l|2004h . The lower and 
upper kinetic temperatures are shown in columns 6 and 7 of 
Tabled As a general rule, the cores at low temperatures are 
well constrained, while those above ~20K are poorly con- 
strained. This is reflected in the larger error bars for poorly 
constrained temperatures, and is not unexpected given the 
insensitivity of NH. 3(1,1) and (2.2) as a temperature probe 
of warmer gas (e.g. iDanbv et al1ll988l ). 

For some sources (e.g. G305. 776-0. 251) the rotational 
te mperature wa s so h igh that the analytic for determined 
bv lTafalla et al.l (|2004f ) was no longer reliable. These sources 
have very small (<8K) lower limits, and no upper limit could 
be derived. 



3.3 Detection Rates 

This Parkes ammonia (1,1) and (2,2) survey targeted a 
total of 2 4 4 sou rces of the 405 in the SIMBA sample of 
iHill et all l|2005h , i.e., 60 per cent. The breakdown of the 



different classes of source targeted and their detection rates 
in both NH 3 (1,1) and (2,2) is given in Table H Of the 244 
sources observed, 138 sources were detected (at a 3-<r thresh- 
old) in NH3(1,1) (56 per cent) including two sources which 



6 T. Hill et al. 
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a Denotes the source class, with MM indicative of MM-only sources, M for methanol mascr associations, R for radio continuum associations and MR for 
sources with both a methanol mascr and radio continuum source. 

^ Indicates the reliability of the baseline, with R— Reliable; BR— Baseline ripple and BL— blended. See section 13.21 

* Denotes sources where the Nti3(2,2) Flux Density is greater than the NH3(1,1). 

* Indicates that the second component of this spectrum (Fig ??) could not be fit. 



had two velocity components in their spectra, and 102 in 
NH3(2,2)(42 per cent). Sources with a methanol maser or 
a radio continuum association (class M and R) have higher 
detection rates (~2/3) than class MM (52 per cent) and MR 
(50 per cent), though this distinction is only slight. Of the 
sources with an NH3(1,1) detection, 74 per cent were also 
detected in the (2,2) transition. Class MR have the highest 
relative detection rate. 



Compared with the number of sources targeted in each 
class, 60 per cent of the sample with active star formation 
(class M, MR and R); i.e., 59 of 96 sources, are detected in 
NH3(1,1) which decreases to 50 per cent for the (2,2) tran- 
sition. In comparison, 50 per cent of the MM-only sources 
are detected in Nfi3(l,l) and only 36 per cent in NHs(2,2). 

Of the remaining 108 sources not reported as detec- 
tions in Nfi3(l,l), 14 sources had obvious detections which 



Table 2. Parameters derived from the NH3(1,1) and NH3(2,2) spectra. Sources are in right ascension order, as per Tablc[T] 



Source 


T rot a 
(K) 


rn lower 
rot 

(K) 


rpuppcr 
rot 

(K) 


T ki „ a 
(K) 


'"plowcr 
kin 

(K) 


rj-,uppcr 

(K) 


(10 13 cm- 2 ) 


N TOT 
NH3 
(10 14 cm- 2 ) 


Rel. 
Grou] 


G269.15-l.13 


28 


18 


63 


42 


22 


- 


5.9 


0.7 


1 


G291. 256-0. 769 


16 


9 


39 


19 


9 


82 


8.6 


1.0 


1 


G291.256-0.743 


22 


19 


28 


29 


22 


11 


8.3 


0.9 


1 


G291. 309-0. 681 


23 


15 


53 


31 


17 


- 


4.8 


0.5 


2 


G291. 576-0. 468 


< 22 


- 


- 


< 22 


- 


- 


4.5 


0.5 


3 


G291.58-0.53 


< 23 


- 


- 


< 23 


- 


- 


4.9 


0.6 


3 


G294.97-1.7 


< 16 


- 


- 


< 16 


- 


- 


4.3 


0.5 


3 


G304. 890+0. 636 


< 9 


- 


- 


< 9 


- 


- 


2.5 


0.7 


3 


G305. 145+0.208 


< 14 


- 


- 


< 14 


- 


- 


4.3 


0.6 


3 


G305. 137+0.069 


< 12 


- 


- 


< 12 


- 


- 


9.7 


1.7 


3 


G305. 192-0.006 


20 


12 


15 


25 


13 


119 


11.7 


1.3 


1 


G305.21+0.21 


37 


26 


66 


71 


37 


- 


12.2 


1.6 


2 


G305. 197+0.007 


14 


10 


22 


16 


11 


28 


15.6 


2.1 


2 


G305.226+0.275 


23 


18 


32 


31 


22 


51 


18.7 


2.1 


2 


G305.228+0.286 


35 


23 


66 


62 


30 


- 


16.6 


2.0 


2 


G305.238+0.261 


18 


12 


31 


21 


13 


18 


10.4 


1.2 


2 


G305.248+0.245 


< 13 


- 


- 


< 13 


- 


- 


8.1 


1.2 


3 


G305.233-0.023 


16 


12 


24 


19 


13 


32 


16.8 


2.1 


2 


G305. 269-0.010 


15 


13 


17 


17 


15 


20 


19.1 


2.5 


1 


G305.355+0.194 


21 


13 


41 


27 


15 


106 


13.2 


1.5 


2 


G305. 37+0.21 


< 16 


- 


- 


< 16 


- 


- 


13.4 


1.7 


3 


G305.362+0.185 


22 


16 


35 


28 


18 


61 


17.8 


2.0 


2 


G305. 361+0.151 


< 26 


- 


- 


< 26 


- 


- 


2.5 


0.3 


3 


G305.538+0.340 


< 14 


- 


- 


< 14 


- 


- 


10.3 


1.1 


3 


G305. 55+0.01 


23 


13 


79 


30 


14 


- 


5.4 


0.6 


2 


G305.552+0.012 


< 26 


- 


- 


< 26 


- 


- 


2.0 


0.2 


3 


G305.561+0.012 


29 


18 


80 


41 


21 


- 


3.2 


0.3 


2 


G305.776-0.251 


17 


2 


- 


19 


2 


- 


1.0 


0.5 


2 


G305.81-0.25 


22 


16 


37 


28 


18 


69 


9.4 


1.0 


2 


G305. 833-0. 196 


12 


7 


26 


13 


7 


37 


6.1 


1.0 


2 


G306.33-0.3 


< 8 


- 


- 


< 8 


- 


- 


6.5 


2.0 


3 


G306.343-0.302 


< 19 


- 


- 


< 19 


- 


- 


8.4 


1.0 


3 


G309.917+0.494 


< 14 


- 


- 


< 14 


- 


- 


8.6 


1.2 


3 


G309.92+0.4 


< 24 


- 


- 


< 24 


- 


- 


3.7 


0.4 


3 


G318.92-0.68 


24 


16 


43 


33 


19 


104 


13.0 


1.1 


2 


G323.74-0.3 


30 


16 


141 


16 


18 


- 


8.0 


0.9 


2 


G332.646-0.647A 


17 


14 


20 


19 


16 


24 


6.3 


0.7 


2 


G332.646-0.647B 


20 


17 


25 


25 


19 


35 


3.1 


0.3 


2 


G332. 695-0.609 


< 10 


- 


- 


< 10 


- 


- 


15.6 


3.1 


3 


G332. 725-0.62 


< 10 


- 


- 


< 10 


- 


- 


9.2 


2.1 


3 


G332. 627-0.511 


17 


12 


28 


19 


13 


41 


6.4 


0.8 


2 


G332. 827-0.552 


18 


12 


32 


22 


13 


52 


16.0 


1.9 


1 


GO.331-0.164 


< 7 


- 


- 


< 7 


- 


- 


12.3 


1.9 


3 


GO.310-0.170 


< 10 


- 


- 


< 10 


- 


- 


11.7 


2.6 


3 


GO.32-0.20 


23 


15 


38 


29 


18 


74 


11.0 


1.3 


2 


Gl. 105-0.098 


32 


27 


40 


53 


39 


83 


5.8 


0.7 


1 


Gl. 13-0. 11 


30 


28 


34 


17 


10 


58 


6.3 


0.7 


1 


G0.549-0.868 


16 


13 


20 


18 


11 


25 


6.5 


0.9 


1 


GO. 627-0. 848 


< 11 


- 


- 


< 11 


- 


- 


13.4 


2.4 


3 


GO.600-0.871 


13 


5 


53 


14 


5 


221 


9.7 


1.4 


2 


G2. 54+0. 20 


12 


9 


19 


13 


9 


23 


10.6 


1.7 


1 


G5.504-0.246 


12 


9 


16 


13 


9 


19 


9.4 


1.6 


2 


G5.90-0.42 


16 


12 


23 


18 


13 


30 


18.9 


2.4 


2 


G5.90-0.44 


23 


16 


40 


30 


18 


84 


10.2 


1.1 


2 


G8.111+0.257 


29 


18 


76 


11 


21 




1.8 


0.2 


2 


G8.127+0.255 


17 


13 


27 


20 


11 


38 


7.2 


0.9 


2 


G8.138+0.246 


17 


13 


26 


20 


14 


36 


8.2 


1.0 


2 


G8.13+0.22 


15 


12 


18 


16 


13 


21 


11.9 


1.6 


1 


G5.948-l.125 


< 9 






< 9 






2.7 


0.7 


3 


G5.962-l.128 


13 


4 




11 


1 




5.2 


0.8 


2 


G5.975-l.146 


21 


15 


35 


26 


17 


61 


6.1 


0.7 


2 


G5.971-l.158 


< 12 






< 12 






5.3 


0.9 


3 


G9.63+0.19 


27 


19 


44 


38 


22 


110 


19.2 


2.2 


2 


G8.68-0.36 


12 


11 


14 


13 


11 


16 


38.9 


6.2 


2 
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Table 2. Cont. 



Source 


T ro t a 
(K) 


i-plowcr 
x rot 

(K) 


rpupper 
1 rot 

(K) 


r-p a 
J^kin 

(K) 


f-p lower 
kin 

(K) 


rpiippcr 
kin 

(K) 


NnH3(1,1) 
(10 13 cm- 2 ) 


tvtTOT 
1N iVif3 
(10 14 cm- 2 ) 


Rel 
Grou] 


PO ftQft n 3fifi 
vjjO.UOU~U.OUU 


1 fi 

1U 


1 3 
10 


9n 

zu 


1 Q 
1 J 


1 ^ 
IO 


24 


3^ f, 

00. u 


4 ^ 

4:.0 


2 


PS fiAA n 3Q^ 
vjj o . u 4: u • o y o 


11 


Q 
O 


1 ^ 

10 


11 


Q 
O 


1 7 

1 1 


Q Q 

y.y 


1 Q 

1 .y 


2 


PS 71 3 n 3fi4 


y 


Q 
O 


1 n 

1U 


Q 

y 


Q 
O 


11 


93 3 
zo.o 


9 
u.z 




ro 7Q-. n 3fi9 

v_jO. iOO-U.OUZ 


12 


Q 


1 7 
1 1 


1 3 
10 


1 n 

1U 


9n 

zu 


1 7 


9 R 
z.o 


2 


r^Q 794 Afll 
vjjO. i Z4:-U.4:U1 


i n 

1U 


7 
1 


1 3 
lO 


1 n 


7 


1 A 
11 


9n 1 
zu. 1 


A « 

4:.U 


1 
1 


PR 70Q n 41 9 
vjjO. luy-u.'iiz 


11 


Q 

y 


14 


12 


Q 

y 


1 fi 

1U 


1^0 

10. y 


9 Q 
z.y 




pr 71 8 n 4in 

VjjO. I 10-U.41U 


i 3 

10 


1 n 

1U 


1 « 
10 


1 A 

14: 


1 1 
1 1 


1 R 
lo 


1 Q R 

iy .0 


3 n 
o.u 


9 
z 


Pin 9R7 n 1 m 

VjjIU. ZO f ~U. 1 1U 


22 


1 ^ 
10 


42 


9Q 
zy 


1 (\ 

1U 


Q3 

yo 


1 R Q 

10. y 


2 1 


2 


Pin 9R4 n 1 9fi 

vjj 1U. ^u4:-U. 1ZU 


9^ 
zo 


1 7 
1 / 


An 

4:U 


33 
OO 


9n 

zu 


RR 

ou 


1 7 4 

1 1 .4: 


9 n 

z.u 


2 


Pin 9RR n 1 97 

VjjIU. ZOO~U. IZ I 


1 

±y 


1 fi 
1U 


9^ 

ZO 


93 
zo 


1 R 
10 


34 

04: 


1 n 7 

1U. f 


1 3 
1.0 




p i n 9Q n 1 4 

VjjIU. Zy-U. 14: 


1 7 
i / 


14 


9 1 

Zl 


9n 

zu 


1 f< 
10 


9fi 
zu 


1 ^ ^ 

10.0 


1 Q 

1 .y 




pin 343 n 1 A9 

v_j 1 U . 04:0 - U . 14:Z 


1 fi 

1U 


11 


97 

Z f 


1 

iy 


12 


3R 
OO 


1 4 3 

14:. O 


1 R 
l.O 


2 


p i n 39 n i ^ 

VjjIU.OZ-U. 10 


24 


1 fi 

1U 


4.0 


33 
OO 


1 R 
10 


120 


1 3 n 
10. u 


1 4 


2 


pin 3^0, n 1 aq A 
vjjiu.ooy-u. 


9fi 
zu 


93 
zo 


9Q 
zy 


37 
O / 


3 1 
01 


44 


3 3 
O.O 


n 3 
u.o 




pin 3^0, n 1 aqr 

vjjiu.ooy-u. i4:y.D 


9^ 
zo 


22 


9Q 
zy 


3^ 
OO 


9R 
zo 


4^ 

4:0 


4 2 


n ^ 
u.o 




p i n «.3 n 33R 

Vjj_LU.UO-U.00-D 


3R 
oo 


9n 

zu 


1 ^7 
10 f 


7^ 
i 


9^ 
zo 




1^1 
10. 1 


1 Q 

1 .y 


2 


pin fi9 n 33 

vjjIU.UZ-U.OO 


1 7 
1 / 


1 3 
10 


24 


9n 

zu 


14 


31 
01 


19 2 


2 4 


2 


pi n <S9 n 3R 

Vjj_LU.U---U.00 


42 


33 


61 


98 


54 




5.5 


0.7 


^ 


PQ RR n 7^ 
vjjy.oo-u. ( o 


14 


12 


1 R 

10 


1 

1U 


12 


21 


1 R R 
lO.O 


9 

z.u 


2 


pi n fi9n n 441 

Vjj _LU. UZU U. Irt-L 


13 


10 


18 


14 


10 


21 


9.3 


1.4 


2 


p 1 1 n7^ n 3R4 

vjjI 1 .U I 0-U.004: 


<v. y 






s y 






12 9 


3 9 
O.Z 


Q 
O 


pi 1 11 n 34 

VJJ _L_L. _L J_ U. 04: 


17 


12 


27 


20 


12 


38 


20.6 


2.5 


2 


pii 117 n 4i ^ 

vjj 11.11 / -U.4:LO 


11 


Q 

y 


1 3 
10 


12 


1 n 

1U 


14 


12 1 


2 2 




P 1 9 RR-Un 4R 

VjjIZ. OOt^U. 4:0 


93 
zo 


1 ^ 
10 


/in 

4:U 


3n 
ou 


1 7 
1 * 


R4 

04: 


14 7 

14:. ( 


1 7 

1 . < 


2 


P 1 9 Q1 A-\-C\ 4Q3 
vjjiz.yi4:-|-u.4:yo 


22 


Q 
O 




9Q 
Z.J 


Q 
O 




9 Q 
z.y 


n 3 
u.o 


2 


pi 1 Q03 n 1 4n 

vjj _L _L . yUO U. 14:U 


12 


9 


16 


12 


9 


18 


21.2 


3.6 




p 1 1 0,3 n i a 

VjjI 1 . yO~U. 14: 


<^ 1 9 
\ iz 






^19 
<^ iz 






13 4 

IO. 4: 


2 1 


Q 
O 


P19 119 n 19^ 

vjj1Z.11Z-U.1Z0 


<. Iz 






^ 19 
<. Iz 






^ 9 

o.z 


n r 

U.o 


Q 
O 


P 1 1 Q49 n 9^fi 
vjj 1 1 . y4:Z-u. zou 


1 Q 

iy 


12 


33 
OO 


93 
zo 


1 3 
10 


OO 


1 r n 
10. u 


2 1 


2 


P19 1R n 19A 

vjjIZ. 10 _ U. LZri 


14 


12 


1 7 
1 / 


1 ^ 

10 


1 3 
10 


1 Q 

iy 


16 4 


9 3 
z.o 




P19 91 n 11Q 
vjj iz . z iu _ u . 1 1 y 


1 7 

1 / 


1 3 

10 


9^ 
zo 


9n 

zu 


14 


3^ 
00 


Q A 

y .4: 


1 1 




p 1 9 43 n n^ 

vjj1Z.4:0-U.UO 


1 ^ 
10 


Q 


34 
1 


1 

ID 


Q 


oy 


^ 7 
0. 1 


n 7 

U. ( 


9 
z 


P 1 9 fiR n 1 R 

Vjj1Z.UO~U. IO 


1 3 
10 


11 


1 

1U 


14 


11 


1 R 
10 


9Q n 
zy .u 


4 4 


2 


P 1 1 QA n R9T3 
vjj 1 1 . y4:~u. uzjj 


11 


1 n 

1U 


12 


12 


11 


1 3 
10 


94 7 

Z4:. ( 


4 ^ 

4:.0 




P 1 1 Q3 n fi1 
vjji i . yo-u.ui 


12 


12 


1 3 

10 


1 3 
10 


12 


1 ^ 
10 


26 1 


4 1 




P 1 9 799 n 91 R 

vjjIZ. ( ZZ-U.ZIO 


\ lO 






<v. IO 






Q R 

y.o 


1 3 
l.O 


Q 
O 


P 1 9 R^ n 99fi 
vjj iz . ooo-u. zzu 


1 fi 
1U 


1 3 
10 


21 


1 Q 

iy 


14 


97 
z / 


1 Q R 
iy .0 


9 ^ 
z.o 


2 


P 1 9 RR 1 ^ n 999 
vjj IZ . ooJ-u. zzz 


1 3 
10 


1 n 

1U 


1 R 
10 


14 


1 n 

1U 


21 


ion 
iy .u 


9 Q 
z.y 


2 


P 1 9 RQ9 n 99fi 
vjj iz . oy z-u. zzu 


11 


Q 

y 


12 


1 1 


1 n 

1U 


1 3 
10 


21 4 


4 1 




p 1 9 on n 9^13 

vjj iz . yu _ u. ZO.D 


1 1 


1 n 

1U 


12 


12 


11 


12 


97 R 
z / .0 


^ 9 
O.Z 




P 1 9 R^Q n 979 
vjj iz . ooy-u. z / z 


22 


14 


42 


9R 

zo 


1 

1U 


Q7 
y / 


16 4 


1 R 
1.0 


2 


P 1 3 R7_l_n 9R 
VjjIO.O i tU.Zo 


40 

4:U 


9 
z 




RR 
OD 


O 

z 




1 r 

l.D 


n 9 
u.z 


9 
z 


p i 9 on n of\ 

vjj iz . y u _ u. zu 


14 


12 


17 


1 

1U 


1 3 
10 


9n 

zu 


3^ ^ 
OO.O 


4 R 

4:.0 


2 


P 1 9 R7R n 99fi 

vjj 1 Z.O 1 O-U.ZiU 


s y 






v. y 






24 4 


fi n 
u.u 


Q 
O 


P 1 9 RQ7 n 9R1 
VjjIZ. oy f ~U. ZOl 


14 


1 n 

1U 


21 


1 fi 

1U 


11 


97 
z / 


28 2 


3 Q 

0. y 


2 


P 1 9 Q1 4 n 9Rn 
vjjiz.y i4:-u. zou 








<. y 






97 3 
z / .0 


4 

0.4: 


Q 
O 


P 1 9 Q3R n 979 
vjj iz . yoo _ u. z ( z 


s y 






^ Q 

<v. y 






1 fi 1 

ID. 1 


3 7 
O. 1 


Q 
O 


P 1 1 4Q 1 4R 

VJJ 1 1 . 4:y~ 1 .4:0 


\ iz 






^ 19 
iz 






7 n 

l .u 


1 1 


Q 
O 


p i a finj~n ni 

vjj14:.UUt^U.U1 


31 
Ol 


1 R 

LO 


yu 


^n 
ou 


21 




1 3 
1U.O 


1 Q 

1 .y 


2 


Pin R4 9 

vjj i u . 04:- z . o y 


22 


12 


Rfi 
ou 


97 
z / 


14 




9 

z.u 


n 3 
u.o 


2 


p i ^Rn n n7Q 
vjjiu.oou-u.u ( y 


/ in 

\ 1U 






^ 1 n 

<v. 1U 






R 9 
O.Z 


1 

1 .u 


Q 




vjjIU.OO-U.UO 


1 A 
14 


1 O 

Iz 


1 Q 

lo 


1 ft 

lu 


1 Q 

lo 


zz 


n q 
y.o 


1 Q 
l.O 


1 

1 


G18.30-0.39 


29 


15 


109 


43 


18 




6.1 


0.7 


2 


G19.61-0.1 


32 


19 


122 


52 


23 




4.9 


0.6 


2 


G16.871-2.154 


14 


13 


17 


16 


14 


19 


28.8 


3.9 


2 


G16.86-2.15 


13 


10 


18 


14 


11 


21 


18.3 


2.7 


1 


G22.36+0.07B 


< 10 






< 10 






10.0 


2.0 


3 


G23.71+0.17 


14 


11 


18 


15 


12 


22 


15.7 


2.2 


2 


G23.43-0.18 


< 7 






< 7 






19.8 


9.0 


3 


G23.268-0.257A 


24 


19 


33 


32 


22 


55 


3.1 


0.3 


3 


G23.409-0.228 


12 


9 


16 


12 


9 


18 


18.6 


3.2 


2 


G23.754+0.095 


< 14 






< 14 






8.8 


1.2 


3 


G29.96-0.02B 


16 


11 


26 


18 


12 


37 


13.4 


1.7 


2 


G29.912-0.045 


14 


12 


18 


16 


13 


21 


22.8 


3.0 


1 
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Table 3. Sources which were observed in NH3(1,1) and (2,2) but 
are non-detections (group 4 sources) or are too confused to fit. 



Table 4. Summary of NH3 detections 



Source Name 


Class 




Source Name 


Class 


1 fe.c 


G269.45-l.47 


MR 


ND 


G305. 519-0. 040 


MM 


ND 


G270.25+0.84 


M 


ND 


G305. 520-0. 020 


MM 


ND 


G284. 271-0. 391 


MM 


ND 


G305. 549+0. 002 


MM 


M 


G284. 295-0. 362 


MM 


ND 


G305. 581+0. 033 


MM 


ND 


G284. 307-0. 376 


MM 


ND 


G305. 605+0. 010 


MM 


M/ND 


G284. 338-0.417 


MM 


ND 


G306. 319-0. 343 


MM 


ND 


G284.35-0.42 


M 


M/ND 


G318. 913-0. 162 


R 


ND 


G284. 345-0. 404 


MM 


ND 


G330. 952-0. 18 


MR 


M 


G284. 341-0. 389 


MM 


M/ND 


G332. 640-0. 586 


MM 


D/M 


G284. 328-0. 365 


MM 


ND 


G332. 701-0. 587 


MM 


ND 


G284. 384-0. 441 


MM 


M/ND 


G332. 777-0. 584 


MM 


M/ND 


G284. 344-0. 366 


MM 


ND 


G332. 794-0. 598 


MM 


ND 


G284. 352-0. 353 


MM 


ND 


GO. 204+0. 051 


MM 


SB 


G287.37+0.65 


M 


ND 


GO. 49+0. 19 


M 


SAT 


G290.40-2.91 


M 


ND 


GO. 266-0. 034 


MM 


SB 


G291.27-0.70 


MR 


D 


GO. 21-0. 00 


MR 


SB 


G291. 288-0. 706 


MM 


ND 


GO. 497+0. 170 


MM 


SAT 


G291. 302-0. 693 


MM 


M 


GO. 240+0.008 


MM 


SB 


G290.37+1.66 


M 


M/ND 


GO. 527+0. 181 


R 


SAT 


G291. 587-0.499 


MM 


ND 


GO. 271+0. 022 


MM 


SB 


G291. 572-0. 450 


MM 


M/ND 


GO. 257+0. 011 


MM 


SB 


G291. 608-0. 532 


MM 


ND 


GO. 83+0. 18 


M 


M/ND 


G291. 597-0. 496 


MM 


ND 


GO. 325-0. 242 


MM 


ND 


G291. 630-0. 545 


MM 


ND 


Gl. 124-0. 065 


MM 


SB 


G291. 614-0. 443 


MM 


ND 


Gl. 134-0. 073 


MM 


SB 


G293. 824-0. 762 


MM 


ND 


Gl. 14-0. 12 


M 


ST 


G293.82-0.74 


MR 


M 


GO. 55-0. 85 


MR 


SB 


G293. 892-0. 782 


MM 


D/ND 


G5. 48-0. 24 


R 


ND 


G293.95-0.8 


MR 


ND 


G6. 53-0. 10 


R 


M/ND 


G293. 942-0. 876 


MM 


ND 


G6. 60-0. 08 


M 


M/ND 


G293. 989-0. 936 


MM 


M/ND 


G6. 62-0. 10 


MM 


ND 


G294.52-1.6 


M 


M/ND 


G5. 97-1. 17 


R 


ND 


G294. 945-1. 737 


MM 


ND 


G10.10+0.72 


R 


ND 


G294. 989-1. 720 


M 


ND 


G9. 966-0. 020 


MM 


M/ND 


G298.26+0.7 


M 


D/ND 


G9. 99-0. 03 


M 


ND 


G299. 02+0.1 


M 


ND 


G10. 001-0. 033 


R 


ND 


G299. 024+0. 130 


MM 


ND 


G10.44-0.01 


M 


D/ND 


G300. 455-0. 190 


MM 


ND 


G9. 924-0. 749 


MM 


ND 


G300. 51-0.1 


M 


ND 


Gil. 948-0. 003 


MM 


ND 


G301.14-0.2 


MR 


ND 


G12.02-0.03 


M 


D/ND 


G302.03-0.06 


MR 


D/ND 


Gil. 902-0. 100 


MM 


M/ND 


G304. 906+0. 574 


MM 


M/ND 


Gil. 861-0. 183 


MM 


ND 


G304. 919+0. 542 


MM 


ND 


Gil. 942-0. 157 


MM 


D/M 


G305. 952+0. 555 


MM 


ND 


G12. 200-0. 003 


MM 


D/M 


G304. 952+0. 522 


MM 


M/ND 


Gil. 956-0. 177 


MM 


ND 


G304. 933+0. 546 


MM 


ND 


Gil. 99-0. 27 


M 


M 


G304. 942+0. 550 


MM 


ND 


G19. 607-0. 234 


MR 


ND 


G305. 201+0. 241 


MM 


M 


G19.70-0.27A 


M 


ND 


G305. 202+0. 230 


MM 


D/M 


G16. 883-2. 188 


MM 


D/M 


G305.20+0.02 


R 


ND 


G21.87+0.01 


MR 


ND 


G305. 200+0. 02 


M 


ND 


G24.450+0.489 


MM 


ND 


G305. 242+0. 225 


MM 


ND 


G23.420-0.235 


MM 


M 


G305. 513+0. 333 


MM 


D/ND 


G23. 319-0. 298 


MM 


ST 


G305. 533+0. 360 


MM 


ND 


G29. 853-0. 062 


MM 


ND 



Denotes the source class. Sec Scctii 



Code indicates whether the source is a: ND — non-detection; M — may be 
a detection, these sources arc typically weak and/or the spectrum is too 
noisy to positively confirm a detection; D — detection too weak compared 
with the noise; S — Strong detection in which cither (B) the hypcrfincs arc 
blended and can't be fit or (A) has absorption features or (^) the 
spectrum is confused with multiple components and a fit will not converge. 
c Two codes appear in this column only if the NH3{2,2) transition is not 
in agreement with the NH3(1,1), with the former following the latter and 
separated by a '/'. 



could not be fit (see Table [3J. These tended to be sources 
quite close to the Galactic Centre and whose spectra dis- 
played very broad lines with a combination of hyperfine 
blending, multiple lines and/or absorption features. A fur- 
ther 11 sources had detections that were too weak (< 3-cr) 
in NHs(l,l) and a further 24 sources are possible detec- 
tions though the noisy spectrum and very low signal-to-noise 



Class 


Sources 


Good Fits 


% of NH 3 (1,1) 




Targeted 


(1,1) 


(2,2) 


with (2,2) 


Total 


244 


138 (56%) 


102 (42%) 


74% 


MM 


148 


79 (53%) 


54 (36%) 


68% 


M 


54 


35 (65%) 


27 (50%) 


75% 


MR 


24 


12 (50%) 


12 (50%) 


100% 


R 


18 


12 (67%) 


9 (50%) 


75% 


M+MR+R 


96 


59 (61%) 


48 (50%) 


81% 



makes is difficult to determine whether there is a detection 
or not. 

Three of the sources targeted have multiple velocity 
components in their spectra: G10.359-0.149; G23. 268-0. 257 
and G332.646-0.647 which are visible in both the NH 3 (1,1) 
and (2,2) transition, though G23.268-0.257B did not pass 
our cuts described in section 13.21 For nine sources, the 
strength of the (2,2) transition is stronger than that of the 
(1,1), indicating that these sources are likely the hottest in 
the sample. It is not possible to determine an accurate col- 
umn density or temperature for these sources, which are an- 
notated by a f in TableQ] These sources do not represent one 
particular class of source. As these transitions were observed 
simultaneously we can rule-out weather effects causing this. 

4 RESULTS 

In this section the results pertaining to the parameters ob- 
tained from fitting the NH3(1,1) and (2,2) spectra (i.e., AV, 
t, Vlsr, Flux density) are discussed, as well as the param- 
eters derived from the fits: Tki n , NH3(1,1) column density 
( N NH 3 (1 1)) and total NH3 commn density (Nj^,). Ta- 
ble [5] presents the mean and median values of each of these 
parameters. 



4.1 Linewidth AV 

The linewidth of the sample ranges from 1.1 to 7.1 kms - 
for the NH3(1,1) transition and 1.3 to 9.2 kms~ for the 
NH 3 (2,2) transition. The mean linewidth of the NH 3 (1,1) 
data is 2.9 kms -1 , with a standard deviation of 1.2kms -1 , 
whilst the (2,2) transition displays a slightly broader 
mean linewidth of 3.1 kms -1 , with a standard deviation of 
1.7 kms -1 . 

IPillai et alj (|2006j ) found linewidths between 0.8 and 
3 kms -1 for their sample of nine infrared dark clouds 
(IRDCs ), whilst for their sampl e of methanol maser selected 
sources. lLongmore et al.l (|2007| ) find linewidths between 0.7 
and 4.6 kms -1 . Both the linewidths of these authors and our 
sample, which is comprised of a cross-section of star-forming 
sources, d isplay linewidths t hat are greater than those re- 
ported by Jiiina et a l. (1999) for their sample of low mass 
cores. IPillai et al.l (|2006l ) attribute the larger linewidths 
of their sample t o vel ocity dispersions and turbulence. 
IChurchwell et al.l (|l990t ) find an average NH3 linewidth 
of 3.1 kms -1 for thei r sample of UC Hn regions, whilst 
ISridharan et al l (|2005l ) find a median linewidth of 1.5 kms 
for their sample of high mass st arless cores and 1.9 km s -1 
for high mass protostellar objects. ISridharan et al.l l|2005l ) in- 
terpret the linewidth difference between the different types 
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Figure 2. Cumulative distribution of the NH3(1,1) (top) and 
NH3(2,2) (middle) linewidths of the sample. Key is as on plot. 
There is a slight offset between the (1,1) and (2,2) transitions. 



of source as an indication of evolutionary status, with more 
quiescent and less evolved cores having smaller linewidths. 

The cumulative distribution plots of the NH3(1,1) and 
NH3(2,2) linewidth (see Fig. [2}, with respect to the class 
of source in the sample, confirms that the MM-only sources 
have the narrowest linewidths of the sample. Class R sources 
are confined to a small region of the linewidth parameter 
space which may simply reflect the small number of sources 
in this sample (see Table [4}. For both the NHs(l,l) and 
(2,2) distributions, the maser sample traces that of the sam- 
ple with active star formation (class M+MR+R) indicating 
that masers are ubiquitous tracers of massive star forma- 
tion. Class MR sources have the broadest linewidths of the 
sample, given that they have two indicators of active star 
formation they could potentially be warmer and/or more 
evolved. 

To test the hypothesis that two, or more, classes 



of source are drawn from the sample parent population, 
Kolmogorov-Smirnov (KS) tests were performed. For both 
the (1,1) and (2,2) transition, the KS-tests for the individual 
classes were inconclusive which is likely due to the low num- 
ber of sources in the individual samples of active sources (M, 
MR and R). When comparing the MM-only sample (class 
MM) with the star formation sample (class M+MR+R), the 
statistics become significant enough to perform robust KS 
tests. The null hypothesis, that the samples are drawn from 
the same population, can clearly be rejected when compar- 
ing the MM-only sample (class MM) with the star formation 
sample (class M+MR+R) - see full red and dotted black 
line in Fig. [2] This results holds for both the (1,1) and (2,2) 
transition, with probabilities smaller than 10 -3 and 7x 10 -5 , 
respectively. 

As the Nfi3(l,l) and (2,2) transition are expected to 
be emitted from the same gas, we would then expect that 
they display similar linewidths. The top two panels of Fig- 
ure [2] suggest that there is a slight increase in linewidth with 
an increase in the NH3 transition. Note that the NH3(1,1) 
plot includes all sources, whilst the NHs(2,2) plot includes 
sources with a Rel. Group equal to 1 or 2, as Group 3 sources 
have only uppe r limits for this tra nsition. This result is in 
agreement with Pillai et al.l (|2006h who also found that for 
some of their IRDCs the (2,2) linewidths were slightly larger 
than the (1,1) linewidths, which they interpreted to mean 
that each transition was not exactly tracing the same gas. 
Broader NH 3 (2,2) linewidths, with respect to the NH 3 (1,1) 
linewidth, could be interpreted as internal heating, but fur- 
ther (mapping) observations of these sources are necessary 
to examine this. 

The bottom panel of Figure [2] is a plot of the difference 
between the NHs(l,l) and (2,2) linewidth. This panel places 
the top and middle panel in context, indicating that approx- 
imately half of the sources have broader (2,2) linewidths and 
the other half have less broad line widths compared with the 
NH3(1,1) transition. This panel also clearly shows that each 
of the classes of source have comparable (2,2) versus (1,1) 
linewidths. If we refer to the sources for which we have ro- 
bust detections (i.e. a Rel. Group equal to 1 or 2), the (2,2) 
linewidths can vary from -2.0 to 2.5 kms -1 broader than the 
(1,1) transition with a median and mean of 0.0 kms -1 . 

For the 45 per cent of sources with Nti3(l,l) linewidths 
broader than the (2,2) transition, the (2,2) linewidth is 50 to 
99 per cent that of the width of the (1,1) transition, with a 
median of 90 per cent of the (1,1) linewidth. The small differ- 
ence in linewidth ensures that we can use the two transitions 
to determine the temperature of the cores. 

Thermal linewidths range from 0.1 to 0.4 kms -1 with 
a median of 0.2 kms -1 , as calculated from the NH 3 derived 
kinetic temperature. The turbulence is clearly dominating 
for all sources, including the MM-only where the linewidth 
is the smallest. 



4.2 Optical Depth r 

The optical depth of the sample ranges from 0.1 to 3.8 for 
the NH 3 (1,1) transition and from 0.1 to 9.4 for the (2,2) 
transition. We caution however that this lower limit is sim- 
ply the minimum value for optical depth as determined from 
CLASS, and should simply be interpreted as optically-thin 
emission. The median optical depth of the NH3(1,1) is 1.4 
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compared with 0.1 for the (2, 2) transition. These op tical 
depth values are consistent with lLongmore et al.l (2007) who 
found that the majority of their maser sources have optical 
depths between 0.3 and 5. There is no obvious difference in 
the optical depth between the different classes of source for 
either of the NH3 transitions. Only 22 per cent of sources 
observed in both (1,1) and (2,2) have greater optical depths 
for the (2,2) transition. 

4.3 V LS R 

The Vlsr of the sample ranges from -87.8 to 113kms~ . 
There are no trends between the source Vlsr and type of 
source in the sample. 

4.4 Flux Density 

The NH3(1,1) flux density of the sample ranges from 0.3 
to 11.60 Jy/beam, both of which are MM-only sources, with 
a median of 2.4 Jy/beam. The different classes of source in 
the sample display a similar range of flux. Class M, MM and 
MR have similar median flux density values of 2.5 Jy/beam, 
whilst class R have a median flux of 1.6 Jy/beam, though the 
sample of active star formation sources (class M, MR and 
R) has a median value of 2.4 Jy/beam. For nine sources, the 
flux density of the (2,2) transition is greater than that of the 
(1,1) - see section [3~3l 

For the NH3(2,2) transition the flux density of the 
sample ranges from 0.3 to 6.4 Jy/beam, with a median of 
1.2 Jy/beam. Each of the four classes of source display simi- 
lar flux ranges for this transition. The median NHs(2,2) flux 
density of each class suggests that the MM-only sources are 
not as bright as the other classes of source. We caution how- 
ever that this difference is small. The NH_3(2,2) flux density 
of the MM-only cores is 75 per cent that of the sample with 
active star formation. 

Cumulative distributions and Komolgorov-Smirnov 
(KS) tests indicate that there is little difference between 
each of the classes in terms of their flux density. 

The peak line intensity of the NH_3(2,2) transition for 
our sources is on average 40 per cent that of the (1,1) tran- 
sition (excluding the nine sources which have greater (2,2) 
flux densiti e s). T his ratio is in excellent agreement with 
IPillai et al.l (|2006l ) who also find that for their sample of 
IRDCs, the peak intensity of the (2,2) is 40 per cent that of 
the (1,1). 

4.5 Kinetic Temperature Tkm 

The kinetic temperature of the sample ranges from 9.0 to 
98.0 K, with a median of 20 K. Despite such a large temper- 
ature range, only 22 per cent of sources have temperatures 
greater than 30 K, whilst only 6 per cent have temperatures 
in excess of 50 K. 78 per cent of sources are cooler than 
30 K. As discussed in section 13.21 and 15.31 the lower transi- 
tion (1,1) and (2,2) NH3 data are usef ul as a thermomete r 
when the temperature is less than 30 K (|Danbv et al.lfl988l ). 
When the temperature is greater than this, these transi- 
tions can not suppo rt accurate temperature determinations. 
iTafalla et ail (|2004h are more conservative and suggest that 
NH3(1,1) and (2,2) actually poorly constrain temperatures 




Figure 3. Cumulative distribution of the column density of the 
different classes in the sample. Key is on the plot. 

greater than 20K. However, our method of modelling, i.e., 
Bayesian inference, factors in the larger error bars for ill con- 
strained temperatures, when deriving masses from SEDs - 
see section [5.2.21 

The MM-only sources are cooler on average than each 
of the individual classes comprising maser and/or radio con- 
tinuum associations as well as the sample of active star for- 
mation (class M+MR+R) - see Table \S\ The cumulative 
distribution plot of the temperature is discussed in detail in 
section 15.11 

IPillai et al.l l|2006l ) found that the te mperature of their 
IRDC sample range from 11-17 K, whilst IChurchwell et al.l 
( 1990) find approximately half of their UC Hn sample to 
range between 15 and 25 K and the other half to have tem- 
peratures in excess of 25 K. 

4.6 Total Column Density 

The total column density of the sources, derived from 
both the NH3(1,1) and (2,2) transitions, ranges from 0.3 
to 15.7 xlO 14 cm~ 2 , with a median column density of 
2.8 xlO 14 cm~ 2 . Figure [3] presents the cumulative distribu- 
tion plot of the column density. The radio sample (class R) 
has the lowest column density on average compared with the 
other classes - see Table [5] 

4.7 Correlations 

Cross-correlating each of the parameters discussed in the 
aforementioned sections reveals very few trends and thus 
relationships between them. There is a slight correlation be- 
tween the flux and the optical depth of a sources, with an 
increase in flux translating to an increase in the optical depth 

of the source. 

Both lLongmore et all (|2008l ) and IPillai et all (|2006l ) 
found that the NH3(1,1) linewidth and kinetic temperature 
of a sources are correlated, with an increase in one leading 
to an increase in the other. lLongmore et al.l (|2008T ) found 
that when they compared these two parameters, the NH3- 
only cores of their sample - those without methanol maser 
associations - were confined to a small and distinct region 
in the plot, which was clearly separate from the parameter 
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space occupied by the methanol maser sources. The NH3- 
only sources were cooler than the maser sources, yet of com- 
parable linewidth. 

As mentioned in section 14.11 iPillai et al.l l|2006l ) found 
that IRDCs in their sample have gre ater linewidt h s tha n 
the low-mass stars in the sample of IJiiina et al.l l[l999). 
Additionally, they fin d that their IRDC s are colder rela- 
tive to th e sample of Jiiina et all (|l999l ) and the UC Hn 
sample of IChurchwell et all (|l990h . A comparison of their 
IRDCs with more evolved examples of massive star forma- 
tion __(e2g;_JJ^_Frn regions and high mass protostellar objects 
e.g. lBeuther et alj |2002) reveals the IRDCs to be cooler than 
these sources and to have narrower linewidths on average. 

We do not see this correlation in our own data, i.e., 
the MM-only sources are not distinct from the other star- 
formation classes with respect to their temperature. Instead, 
we find that the MM-only sample has comparable temper- 
atures to sources with a methanol maser and/or radio con- 
tinuum sources, coupled with smaller linewidths. The eight 
dark clouds for which we have NH3(1,1) and (2,2) obser- 
vations are not c onfined to cooler te mperatures or smaller 
linewidths as per IPillai et al. I (|2006h . nor is there a trend 
between their NH3(1,1) linewidth and Tki n - The lack of 
correlation between the NH3(1,1) and Tkin of our sam- 
ple may be attributed to two factors. Firstly, our data are 
likely subject to beam dilution with low-resolut ion (58 arc- 
sec), c ompared with the 40 arcsec resolution of IPillai et al.l 
(2006) and the in t erfero metric observations (< 11 arcsec) 
of Longmore et al.l {2008). Additionally the combination of 
NH3(1,1) and (2,2) provides reliable temperatures up to 
~ 30 K, so we ar e insensitive to t he gr eater temperature 
range coverage of lLongmore et aL (2008) who used higher 
transition a mmonia data to d eterm ine their temperatures. 

While lLongmore et all <|2006l ) were able to probe 
warmer regions better than us, due to higher transition NH3 
data which they also measured, when it comes to the cold 
sources of particular interest here, this is not an issue. Both 
studies are sensitive to probing the temperature in the cold- 
est gas, where only the (1,1) and (2,2) lines are significantly 
excited. 



5 ANALYSIS 

5.1 Previous SED modelling 

iHill et all (|2009h performed spectral energy distribution 
modelling of ~ 180 sources of the SIMBA sample (see sec- 
tion [Tj), using the Bayesian inference method of fitting (e.g. 
iPinte et al.ll2008T ). This method of fitting considers the po- 
tential correlations between parameters to produce quanti- 
tative estimates of the range of validity of key parameters 
(temperature, mass, luminosity) extracted from SED fitting. 
As th e modelling proc edure is outlined in detail in section 
3.1 of iHill et alj ( 20091 ) we simply summarise the procedure 
here. 

The sources were modelled according to a two- 
component model which denotes a central warm core sur- 
rounded by a cold dust envelope. The hot component of 
this model is assumed to radiate as a blackbody whilst the 
cold component account s for optically t hin emission from 
the dust (see Equation l. lHill et alj|2009h . The sources were 



then fit for four free parameters: i?hot, Thot, Tcoid and M co id. 
Those sources without mid-infrared emission (i.e, a hot com- 
ponent) were fit for the cold component only and thus two 

free parameters: T co id and M co id. 

It was clear from SED modelling (IHill et al.ll2009l ) that 
an absence of far-infrared data, where the peak of the dust 
emission lies, hinders accurate determinations of the source 
temperature. Additionally as the mass and luminosity of a 
source are highly correlated with that of the temperature, 
assessing the evolutionary status of a source from SED fit- 
ting alone is difficult. Greater observational constraints were 
necessary in order to facilitate further SED fitting and anal- 
ysis. 



5.2 SED modelling revisited 

5.2.1 Comparison of temperature derivation methods 

These ammonia data provide an independent, and more ac- 
curate, determination of the source temperature. Following 
the assumption that the kinetic temperature is equivalent 
to the dust temperature, it is the n poss ible to revisit our 
previous SED modelling. iLi et al.l (|2003l ) showed that the 
gas temperature of their sources i.e., the kinetic tempera- 
ture as derived from NH3 observations, was within a few 
K of the dust temperature. According to their observations, 
Tkin and Td us t are expe cted to be similar in cold regions. 
ISchnee fc Sargentl (|2007l ) also find excellent agreement be- 
tween the dust and gas tempera ture of their st a r less core 
in Taurus. Cont r ary to this, both lMolinari et al.l (1 19961 ) and 
ISridharan et al.l l|2002h find discrepancies when comparing 
dust temperatures and kinetic temperatures derived from 
NH3. However, both these authors derived their dust tem- 
peratures from IRAS data which is subject to poor reso- 
lution and the emission is likely optically thick. When the 
dust is highly optically thick large dif ferences between the 
gas a n d dust tempera t ures c an occur |Kruegel fc Walmslevl 
ll984l ). lMolinari et al.1 (Il996l ) themselves caution that IRAS 
fluxes alone are insufficient for proper estimates of the dust 
temperature. 

Of the ~ 180 sou rces that were originally SED mod- 
elled (|Hill et al.ll2009T) . 82 were also observed with ammo- 
nia, though 30 of these sources only have upper limits to 
their temperature. Figure [4] compares the kinetic tempera- 
ture (Tkin) of each source, as derived from the NH3 obser- 
vations, with the dust temperature as estimated from our 
previous SED modelling (Tsed). This figure contains the 
52 source with reliable continuum and ammonia observa- 
tions for estimating both temperatures (Rel. Group 1 and 
2, in Table [2}, thus allowing a statistical comparison of both 
methods. The red line on the plot indicates where Tsed and 
Tkin are equal. 

Figure [4] shows that there are no systematic biases in- 
troduced in either of the methods used to determine the 
temperature estimates. If SED fitting was biasing the re- 
sultant temperature of a core, then these data would all be 
above, or below, the red line. The SED method, although 
not strongly constraining the temperature, appears a ro- 
bust method for obtaining a first order unbiased estimate of 
the core temperature. In this section, we perform the SED 
modelling again, this time using the kinetic temperature as 
derived from NH3. 
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Table 5. Mean and median as derived from the ammonia spectra for each of the source classes. For a breakdown of the number of 
sources within each class, refer to Table l4l 
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Figure 4. Comparison of the temperature as derived from SED 
modelling with that derived from the NH3 observations. The red 
line indicates where these two parameters would be equal. The 
error bars are indicated by the blue lines. 



5.2.2 SED modelling with NH 3 

SED modelling was performed using the Bayesian inference 
method, following the procedure summ arised above (sec- 
tion [5TTJ and outlined in lHill et al.l (|2009l . section 3). Rather 
than fitting for the temperature (T co id), we fixed the dust 
temperature to the kinetic temperature, as derived from 
NH3, and fit solely for the mass of the source (as well as the 
hot parameters, for those sources with mid- infrared data), 
assuming t he sa me dust properties used previously. As per 
iHill et all (|2009h we have assumed the near distance for all 
sources with a distance ambiguity. 

The probability distributions of the temperature were 
set to Gaussian distributions defined by the value of, and 
uncertainties with, the NH3 temperature (Tki n ) - see Ta- 
ble Q] For those sources where only an upper limit value to 



the temperature has been determined (Rel. Group 3) we as- 
sume a uniform probability distribution between 2.73 K and 
the upper limit. The estimated range of validity for the mass 
and the luminosity derived from this method then not only 
accounts for the uncertainty on the temperature but also 
considers the correlations between parameters. 

Figure [5] is a plot of the probability curves for the tem- 
perature and m ass of a source, c omparing our previous SED 
fitting method (|Hill et al.ll2009l ). with that of the SED fit- 
ting performed in this paper using NH3 derived tempera- 
tures. Figure [5] illustrates how kinetic temperatures allow 
tighter constraints on the mass of a source, as is expected. 
This source, G10. 288-0. 127 a MM-only source with no mid- 
infrared association which was consequently fit with a single 
cold component , corresponds to the bottom row of Fig. 1 of 
IHill et all (120091 ). This source was chosen to illustrate how 
a source with limited sampling of the SED can be well con- 
strained by known or constrained temperatures. It is clear 
that the NH3 data (red curve) provides greater constraints 
on the temperature (left-hand plot), especially for those 
sources with few data points, compared with deriving the 
temperature from SED fitting (black curve). The range of 
validity for the temperature is now constrained to a smaller 
portion of the parameter space. It is clear that accurate tem- 
peratures are necessary to derive accurate sources masses. 



5.3 SED results 

Table [6] presents the range of validity for the tempera- 
ture,mass and luminosity of each of the 52 sources which 
were re-modelled using their NH3 temperature in the SED 
fits. The cumulative distribution plots of the temperature, 
mass and luminosity are presented in Figure [5] Note that 
all sources, regardless of their reliability flag (Table [5]) were 
included in this Figure. 

An interesting result evident from the probability distri- 
butions of all of the sources re-modelled is the complemen- 
tarity of the SED fitting and NH3 observations in terms of 
parameter determinations. Even when the uncertainty on 
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Table 6. Parameters resulting from SED modelling of 52 SIMBA sources with NH3 temperatures reported in Table [2] For each of the 
temperature, mass and luminosity, the range corresponding to a 1-sigma (68 per cent) probability of occurrence is presented, with the 
m i„ and max values representing the lower and upper values of this range, respectively. 



Ident Fit Temperature Mass Luminosity 

Source Name tracer Type T coldmiii T coldmax M min M max L m in Lmax 

K K M M Q L L q 



G305. 833-0.196 


mm 


SINGLE 


26 


28 


7 


6E+01 


1 


.0E-i 


-02 


3 


.7E+03 


1 


.6E+04 


G323. 74-0.3 


m 


SED 


15 


37 


2 


.OE+02 


2 


,8E^ 


-03 


5 


.3E+03 


1 


.OE+05 


GO.32-0.20 


mr 


SINGLE" 


17 


31 


3 


.3E+03 


8. 


.7E+03 


3 


.4E+04 


6 


.7E+05 


Gl. 105-0.098 


mm 


SINGLE" 


28 


47 


5 


.3E+02 


1 


,4E^ 


-03 


1 


.OE+05 


9 


.8E+05 


Gl. 13-0.11 


r 


SED" 


34 


47 


3 


.1E+02 


3 


.8EH 


-03 


4 


.6E+05 


1 


.4E+06 


G0.549-0.868 


mm 


SINGLE" 


15 


24 


2 


.5E+01 


5. 


0E-I 


-01 


6 


.1E+01 


5 


.3E+03 


GO.600-0.871 


mm 


SINGLE" 


11 


42 


1 


.9E+01 


1 


.3EH 


-02 


8 


.9E+01 


7 


.2E+04 


G2. 54+0.20 


m 


SINGLE 


17 


26 


2 


.OE+02 


4. 


6EH 


-02 


5 


.7E+02 


2 


.4E+04 


G5.504-0.246 


mm 


SINGLE" 


10 


19 


1 


.9E+03 


5. 


0EH 


-03 


8 


.3E+02 


3 


.4E+04 


G5.90-0.42 


rn 


SED" 


15 


29 


5 


.3E+02 


1. 


9E-I 


-03 


5 


.3E+03 


7 


.2E+04 


G8.111+0.257 


mm 


SINGLE" 


21 


83 


6 


.1E+00 


3 


.3EH 


-01 


5 


.7E+02 


6 


.7E+05 


G8.127+0.255 


mm 


SINGLE" 


14 


32 


5 


.7E+01 


2. 


0EH 


-02 


2 


.7E+02 


5 


.OE+04 


G8.138+0.246 


mm 


SINGLE" 


15 


31 


1 


.1E+02 


4. 


0EH 


-02 


5 


.7E+02 


7 


.2E+04 


G8.13+0.22 


mr 


SED" 


14 


20 


9 


.3E+02 


2 


.5EH 


-03 


3 


.7E+03 


3 


.4E+04 


G5.962-l.128 


mm 


SINGLE" 


8 


31 


8 


.1E+00 


7. 


6E^ 


-01 


6 


.6E+00 


1 


.6E+04 


G5.975-l.146 


mm 


SINGLE" 


17 


50 


7 


.1E+00 


2 


.5E+01 


1 


.9E+02 


7 


.2E+04 


G9.63+0.19 


mr 


SED" 


16 


36 


2 


.OE+02 


1 


.1EH 


-03 


2 


.5E+03 


5 


.OE+04 


G8.68-0.36 


mr 


SINGLE" 


11 


16 


6 


6E+03 


1 


.1EH 


-04 


3 


.7E+03 


3 


.4E+04 


G8.686-0.366 


m 


SINGLE" 


16 


23 


1 


.1E+03 


2 


,2E^ 


-03 


5 


.3E+03 


7 


.2E+04 


G10.287-0.110 


mm 


SINGLE" 


12 


35 


5 


.7E+01 


3 


1E+02 


2 


.7E+02 


7 


.2E+04 


G10.284-0.126 


m 


SED" 


18 


45 


4 


3E+01 


2. 


3EH 


-02 


1 


.2E+03 


3 


.4E+04 


G10.288-0.127 


mm 


SINGLE" 


19 


32 


2 


.8E+01 


6. 


6EH 


-01 


3 


.9E+02 


1 


.6E+04 


G10.29-0.14 


mr 


SED" 


16 


25 


2 


.7E+02 


7 


.1EH 


-02 


3 


.7E+03 


2 


.4E+04 


G10.343-0.142 


m 


SINGLE" 


15 


35 


5 


.0E+01 


1. 


.5EH 


-02 


2 


.7E+02 


5 


.OE+04 


G10.63-0.33B 


mm 


SINGLE" 


16 


62 


1 


.7E+02 


1 


.2EH 


-03 


1 


.1E+04 


4 


.3E+06 


G10.62-0.33 


m 


SED" 


14 


28 


9 


.3E+02 


3 


.8EH 


-03 


5 


3E+03 


7 


.2E+04 


G9.88-0.75 


r 


SED" 


13 


20 


1 


.2E+03 


3 


.3EH 


-03 


2 


.5E+03 


1 


.6E+04 


G10.62-0.38 


mr 


SINGLE" 


12 


32 


6 


6E+03 


3 


.5EH 


-01 


2 


.4E+04 


2 


.1E+06 


Gil. 11-0.34 


r 


SED" 


13 


28 


8 


.1E+02 


3 


.8EH 


-03 


3 


.7E+03 


7 


.2E+04 


Gil. 117-0.413 


mm 


SINGLE" 


11 


14 


4 


.6E+02 


8. 


TEH 


-02 


8 


.9E+01 


3 


.7E+03 


G12.88+0.48 


m 


SINGLE" 


14 


28 


9 


.3E+02 


3 


.3EH 


-03 


1 


.7E+03 


1 


.OE+05 


G12.914+0.493 


mm 


SINGLE" 


9 


36 


5 


.7E+01 


4. 


.6EH 


-02 


8 


.9E+01 


7 


.2E+04 


Gil. 903-0.140 


mr 


SINGLE" 


9 


17 


5 


.3E+02 


1 


6E-I 


-03 


1 


.3E+02 


1 


.1E+04 


G12.18-0.12A 


m 


SINGLE" 


14 


19 


2 


.5E+03 


4. 


.3EH 


-03 


3 


.7E+03 


3 


.4E+04 


G12. 216-0. 119 


mm 


SINGLE 


16 


31 


2 


.2E+03 


5. 


.7EH 


-03 


2 


.4E+04 


6 


.7E+05 


G12.43-0.05 


r 


SED 


16 


30 


2 


.2E+03 


8. 


.7EH 


-03 


1 


.1E+04 


3 


.2E+05 


G12.68-0.18 


m 


SED" 


11 


17 


1 


.9E+03 


4. 


.3EH 


-03 


2 


.5E+03 


1 


.1E+04 


G11.94-0.62B 


mm 


SINGLE" 


11 


13 


1 


.6E+03 


2. 


.5EH 


-03 


5 


.7E+02 


5 


.3E+03 


Gil. 93-0.61 


mr 


SED" 


13 


15 


1 


.6E+03 


3 


3E+03 


1 


.7E+03 


7 


.7E+03 


G12.90-0.25B 


mm 


SINGLE" 


11 


12 


7 


.1E+02 


1. 


.1EH 


-03 


1 


.9E+02 


3 


.7E+03 


G13.87+0.28 


m 


SED 


10 


39 


i 


.6E+02 


6. 


,6E^ 


-03 


7 


.7E+03 


3 


.2E+05 


G12.859-0.272 


mm 


SED" 


15 


40 


1 


.3E+02 


8. 


TEH 


-02 


1 


.2E+03 


3 


.4E+04 


G12.90-0.26 


m 


SED" 


14 


19 


1 


6E+03 


3 


.8EH 


-03 


5 


3E+03 


2 


.4E+04 


G14.60+0.01 


mr 


SED" 


11 


31 


1 


.3E+02 


9. 


.3EH 


-02 


3 


.9E+02 


1 


.1E+04 


G10.84-2.59 


r 


SINGLE" 


12 


30 


1 


.3E+02 


6. 


TEH 


-02 


1 


.9E+02 


3 


.4E+04 


G16.58-0.05 


m 


SED" 


14 


21 


8 


.1E+02 


2. 


.2EH 


-03 


2 


.5E+03 


1 


.6E+04 


G18.30-0.39 


r 


SED 


14 


36 


1 


.7E+02 


1 


.4EH 


-03 


5 


.3E+03 


7 


.2E+04 


G19.61-0.1 


m 


SED" 


13 


36 


1 


.3E+02 


1. 


TEH 


-03 


5 


.7E+02 


5 


.OE+04 


G16.86-2.15 


m 


SED" 


11 


17 


8 


.1E+02 


2 


.2EH 


-03 


8 


.3E+02 


3 


.7E+03 


G23.71+0.17 


r 


SED" 


12 


20 


1 


.6E+03 


5. 


0E-I 


-03 


7 


.7E+03 


7 


.2E+04 


G29.96-0.02B 


mr 


SED" 


13 


31 


2 


.2E+03 


1 


0E-I 


-04 


1 


.1E+04 


4 


.6E+05 


G29.912-0.045 


mm 


SINGLE" 


11 


20 


1 


.9E+03 


3 


.8EH 


-03 


3 


.7E+03 


5 


.OE+04 



denotes those sources that were fit with IRAS upper limits - see iHill et alj feOOSft . 
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T (K) M (M sun ) 

Figure 5. The probability curves of one source (G10. 288— 0.127, corresponding to the bottom row of Fig. 1 Iffill et al.ll2009t) comparing 
the two SED methods. The black line indicates the probability curve corresponding to a SED fit temperature l|Hill et alj|2009h whilst 
the red probability curve is indicative of the probability when the temperature is fixed using the ammonia derived Tk; n . As can be seen, 
constraints on the temperature translate as tighter probability curves for parameters resulting from SED modelling, i.e., mass. 



the kinetic temperature derived from NH3 is large, the con- 
straints from the SED fitting and the NH3 analysis are very 
complementary. SED modelling provides good constraints 
on the upper limit of the temperature when the ammonia 
data are ill-constrained due to the absence of higher tran- 
sitions (see section [372]) . On the other hand, NH3 provides 
good constraints on the lower limit to the temperature when 
the ambiguities from SED modelling become very large. 

Comparison of the temperature-fit probability curves 
ijHill et aLl llofih with those of the probability curves pro- 
duced here, indicates that, in most instances, the mass prob- 
ability peaks in the same place for both methods. In all in- 
stances, the probability curve of the mass for each individual 
source is more tightly constrained using the NH3 tempera- 
tures, as is demonstrated in Fig. [5] 

The cumulative distribution plot of the temperature 
(Fig. [5J left) suggests that class MM sources are slightly 
cooler than the other classes of source. This difference re- 
flects the fact that 9 of the 11 sources in Rel. Group 3, 
and thus sources with upper limits, are MM-only sources 
which effectively add a 'bump' to the MM-only curve at low 
temperatures. However, the difference between the MM-only 
sources and those with signatures of massive star formation 
is not statistically significant, with a KS test probability of 
0.4. If we do not include the sources in Rel. Group 3, the 
distribution of the temperature shows little distinction be- 
tween the different classes of source in the sample for this 
parameter with a KS probability of 0.98. 

The cumulative distribution plot of the mass (Fig. H3 
middle) indicates that the MM-only sources are the least 
massive of the sample, whilst class MR are the most massive. 
The difference in the mass is more significant than for the 
temperature, though the KS test, with a probability of 0.08, 
does not allow us to rule on the null hypothesis that MM- 
only sources are drawn from the same population as that 
star formation sources. This is likely a result of the relatively 
small number of sources modelled. 



The cumulative distribution plot of the luminosity 
(Fig. [6] right) reveals the MM-only sources to be marginally 
less luminous than the other classes of source, which is not 
surprising given that the mass and luminosity are correlated 
and the temperature of the sources is slightly smaller. 

5.4 Luminosity vs Mass Diagram 

As per iHill et al l l|2009h . we have drawn a luminosity vs. 
mass (hereafter M — L diagram) diag ram for the sources 
which we have re-modelled with SEDs. lAndre et al] (2000) 
proposed that the M — L diagram was a useful diagnostic 
tool for class and class I low-mass protostars, which pro- 
vides insight into a source's evolutionary status. Irlill et all 
(2009) showed that care should be taken with interpretation 
of M — L diagrams drawn from SEDs. 

Due to the low number of sources which were re- 
modelled in each class of source, little distinction can be 
discerned for the different classes of source nor can proposed 
evolutionary tracks be drawn from this diagram. The M — L 
diagram can be found in Figure [7] 

5.5 Determining a mass for NH3 sources not SED 
modelled 

Of the 138 sources deemed to have good NH3(1,1) fits re- 
liable kinetic temperatures were derived for 128 (see sec- 
tion [3T3J) . Of these 128 only 52 NH3 sources were originally 
SED modelled by us. Consequently, only 52 NH3 sources 
have enough information to draw spectral energy distribu- 
tions and determine their mass. 

Assuming that the NH3 temperat ure is equiva- 
lent to that of the dust temperature (|Li et al.l 120031 ; 
ISchnee fc Sargenu 120071 ) . we can estimate the mass of the 
remaining sources which have reliable NH3 and millimetre 
continuum data. Here we ignore sources which are lacking 
millimetre continuum or distance information. 
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Figure 6. Cumulative distributions for the SIMBA sources that were re-modelled with SEDs using their NH3 temperature. The key on 
the left plot indicates the different classes of source in the sample, left: distribution of temperature, middle: distribution of mass, right: 
distribution of luminosity. 




M (Msun) 

Figure 7. Mass-Luminosity diagram for the different classes of 
source in the sample. Each of the points represents the median 
value (i.e. the values at which the cumulative probability distri- 
bution is 0.5) of the mass and luminosity of the individual sources 
(green cross: mm, red plus sign: m, blue triangles:mr and pink cir- 
cles: r). The contour levels represent the regions enclosing 68 per 
cent of the Bayesian probability for a given class of sources (green 
full line: mm, red dashed line: m, blue dot-dash line: mr and pink 
dot-dot-dash line: r). The full and dashed straight lines depict 
the M = L and M = L° 6 relation s, respectively. The same figure 
can be found in Irlill et al. | ||2009| , Fig. 3) for sources which were 
previously SED modelled. 

The 65 NH3 sources for which we derive a mass from 
their NH3 temperature, accordin g to the procedur e outlined 
in section 4.1, and equation 1, of Irlill et alj (120051 ). are pre- 
sented in Table [7] Sources which have reliable temperature 
information and an estimate of the error associated with the 
temperature have both a minimum and maximum estima- 
tion of the mass in columns 3 and 4. Those sources for which 
a mass was derived from an upper limit temperature value 
do not have a temperature range and are lacking values in 
these columns. These masses should be treated with caution 
and are an estimate of the mass only. 

Though the mass drawn from this method is derived 
solely from the millimetre continuum flux, rather than a 



SED sampled at various wavelengths, it provides a useful 
approxima tion to th e sourc e mass in the absence of fur- 
ther data. iHill et all (|2005T ) showed that the contribution 
of free-free emission to the millimetre flux is expected to be 
negligible compared with the dust emission. We recognise 
the potential for contributions from free-free e mission that 
may be missed from the interferometric survey (|Walsh et al.l 
used to draw these conclusions and we caution inter- 
pretation of these results (cf. lLongmore et al. I l2009h . 



6 DISCUSSION 

6.1 Physical properties of the clumps 

Previous work has shown the MM-only cores of our sam- 
ple to be excellent candidates for young massive protostars. 
Their characteristic physical properties, such as tempera- 
ture and mass, suggested that these objects are indicative 
of the earliest stages of massive star formation. These con- 
clusions, which were drawn from (sub) millimetre continuum 
emission observations and spectral energy distribution mod- 
elling, are not without their caveats. SED modelling is sub- 
ject to ambiguities and is heavily reliant on well constrained 
temperatures in order to produce a reliable mass estimate. 
In addition, it is not yet known whether the MM-only core 
is indeed forming massive stars. 

Using the kinetic temperatures as derived from this NH3 
line study, we are now in a position to ascertain which phys- 
ical properties of the MM-only cores are similar to those 
of sources with known star formation activity, and which 
physical properties differ. 

Analysis, in the previous sections, has revealed the MM- 
only class to have smaller NH3 linewidths on average than 
class M, MR and R - those sources with a methanol maser 
and/or radio continuum association. This is true for both 
ammonia transitions. SED modelling, using the Bayesian 
inference method, in conjunction with well constrained ki- 
netic temperatures, as derived from ammonia, reveals the 
MM-only cores to be the least massive of the sample. 

Interestingly, class MM sources display many similar- 
ities in the distributions of the other physical parameters 
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Table 7. Mass of the 65 NH 3 cores that could not be SED mod- 
elled. A minimum and maximum mass are determined from the 
errors on the temperature (see Table [TJ. For those sources in Rcl. 
group 3 an error on the temperature could not be obtained, and 
in this instance there is no corresponding error on the mass given 
in this table. In a few instances, no physically meaningful max- 
imum temperature could be derived, and these sources will not 
have a minimum mass in this table. 



Name 




M 


M min 


Mmax 


G269.15-l.13 


4 


.Oe+02 


1.8e+01 


1 


.0c+03 


G291. 256-0. 769 


4 


.6e+02 


9.0e+01 


1 


.6e+03 


G291. 256-0.743 


5 


,3e+02 


3.4e+02 


7 


.2e+02 


G291. 309-0.681 


6 


.lc+02 


8.1c+01 


1 


.4e+03 


G291. 576-0.468 


1 


le+03 






- 


G291. 58-0.53 


4 


.3e+03 






- 


G294.97-1.7 


1 


,5e+02 






- 


G305.145+0.208 


1 


.3c+02 






- 


G305.137+0.069 


5 


,3e+02 






- 


G305.192-0.006 


2 


.3e+02 


4.5e+01 


6 


.4e+02 


G305.197+0.007 


4 


,6e+02 


2.3e+02 


8 


.6e+02 


G305.21+0.21 


5 


.3e+02 


1.2e+01 


1 


.3e+03 


G305.226+0.275 


1 


.lc+02 


6.7c+01 


1 


.8e+02 


G305.228+0.286 


9 


.3e+00 


2.5e-01 


2 


.3c+01 


G305.238+0.261 


3 


.le+02 


1.2e+02 


6 


.6e+02 


G305.248+0.245 


4 


,6e+02 






- 


G305.233-0.023 


1 


.Oe+02 


5.5e+01 


1 


.8e+02 


G305.269-0.010 


6 


.le+02 


5. lc+02 


7 


.6e+02 


G305.362+0.185 


3 


.5e+02 


1.4e+02 


6 


.2e+02 


G305.361+0.151 


3 


.lc+02 






- 


G305.37+0.21 


9 


.3e+02 






- 


G305.55+0.01 


1 


.lc+02 




3 


.4e+02 


G305.552+0.012 


2 


.Oe+02 






- 


G305.561+0.012 


2 


.7c+02 




6 


.7e+02 


G305.776-0.251 


1 


.lc+02 




5 


.6e+04 


G305.81-0.25 


9 


3e+02 


3.7e+02 


1 


.8e+03 


G306.33-0.3 


1 


.le+02 






- 


G306.343-0.302 


3 


,8e+01 






- 


G309. 917+0. 494 


1 


.3e+02 






- 


G309.92+0.4 


1 


.4e+03 






- 


G318.92-0.68 


1 


.7c+02 


4.9e+01 


3 


.5e+02 


G332.695-0.609 


1 


.9e+03 






- 


G332.725-0.62 


3 


.le+02 






- 


G332.627-0.511 


1 


.Oe+02 


3.8e+01 


1 


.7e+02 


GO. 331-0.164 


2 


,5e+03 






- 


G0.310-0.170 


5 


.3e+02 






- 


GO. 627-0. 848 


5 


.0e+01 








G5. 90-0. 44 


3 


.5e+02 


1.2e+02 


7 


.lc+02 


G5.948-l.125 


6 


,6e+01 








G5.971-l.158 


8 


.7e+01 








G8.644-0.395 


2 


.7e+02 


1.3e+02 


1 


.5e+02 


G8. 713-0.364 


1 


lc+03 


7.8e+02 


1 


.4e+03 


G8.735-0.362 


1 


lc+03 


5.7e+02 


1 


.6e+03 


G8.709-0.412 


i 


.Oe+02 


2.8e+02 


7 


.le+02 


G8.724-0.401 


2 


,3e+02 


1.4e+02 


1 


.6e+02 


G8. 718-0.410 


8 


,7e+01 


6.1c+01 


1 


.3c+02 


G10.32-0.15 


1 


,5e+02 


3.9e+01 


3 


.5e+02 


G10.620-0.441 


3 


.5e+02 


2. lc+02 


6 


.4e+02 


Gil. 075-0.384 


1 


.le+03 








Gil. 93-0.14 


3 


,5e+02 








Gil. 942-0.256 


2 


.Oe+02 


7.0e+01 


4 


.3c+02 


G12. 112-0. 125 


3 


,3e+03 








G12. 722-0.218 


6 


.lc+02 








G12.885-0.222 


1 


.7e+02 


9.5e+01 


2 


.9e+02 


G12.892-0.226 


1 


,5e+02 


1.2e+02 


1 


.8e+02 


G12. 878-0.226 


2 


.Oe+02 








G12. 897-0.281 


8 


.7e+01 


4.6e+01 


1 


.6e+02 



Table 7 - continued 



Name 


M 


M min M max 


G12. 914-0.280 


1.7e+02 




G12.938-0.272 


1.5e+02 




Gil. 49-1.48 


5.3e+02 




G16.580-0.079 


4.0e+02 




G22.36+0.07B 


2.5e+03 




G23.43-0.18 


1.3e+04 




G23.409-0.228 


1.2e+03 


6.6e+02 1.9e+03 


G23.754+0.095 


3. le+02 





compared with the other classes of source. They have sim- 
ilar luminosities, NH3 flux densities hence brightness level, 
and column densities as the other classes. The simultane- 
ous analysis of the NH3 and continuum data shows that the 
distribution of the temperature (Fig. EJ) for the MM-only 
sources is not very distinct from that of the other classes, 
i.e., those with known star formation activity, but seems to 
extends to slightly smaller temperatures. 

6.2 Are they gravitationally bound? 

From SED modelling, it is clear that both NH3 and con- 
tinuum observations provide complementary constraints on 
SED derived source parameters. Even if each of these con- 
straints, taken individually, provide limited information on 
the nature of the cores, the combination of these constraints 
reveal several trends that give insight into the cores and 
their evolutionary status: Class MM are smaller (jHill et al.l 
2005), less massive and more quiescent (smaller turbulent 
linewidths) than cores with signatures of star formation. 

With this information, we can ascertain the gravita- 
tional equilibrium of the various classes of source in the sam- 
ple. A plot presenting the virial mass M v i r i a i = AV 2 R/G, as 
a function of the gas mass derived from continuum observa- 
tions (i.e. SIMBA), M cont muum, is presented in Fig. [8] This 
plot contains all sources which have a reliability flag (Rel. 
Group) of 1 or 2 in Table [2] though as a check, by includ- 
ing the sources with a reliability flag of 3, the correlation 
coefficient and the line of best fit do not differ significantly. 

As can be seen in Fig [8] there is a strong correlation 
between the virial mass and gas mass of a source, with a 
correlation coefficient of 0.74. As the continuum mass of a 
source increases so too does the virial mass. Interestingly, 
the distribution of the virial parameter is very similar for 
the MM-only sources and the active star formation sample 
(class M+MR+R) suggesting a very similar gravitational 
state, despite the lower masses of the MM-only sources. 

The virial parame ter a = M vil .iai/A^contmuum 
|Bertoldi fc McKeel Il992l ) has a median value of 0.95 
(see Fig. [8] right panel), suggesting that the clumps are 
virialised. The precise value must be interpreted with 
caution given the uncertainties on both the virial mass (in 
particular on the size of the clumps) and on the continuum 
mass (e.g. uncertain dust emissivities), but it is clear that 
the clumps, on average, do not differ significantly from 
being virialised. 

There is however a distribution of the virial parameter 
as a function of the source mass, ranging from ~ 0.1 to ~ 10. 
The virial mass appears to vary as M v i r iai oc A/ con ti nuum , in- 




Continuum mass [M SU n] Virial parameter 

Figure 8. Left: The virial mass of the source as a function of the continuum mass. The blue line is where the two masses are equal, 
whilst the red line is the fit to the data, which has a correlation co-efficient of 0.74. The average error bars are as indicated on the plot, 
assuming a conservative 50 per cent error on the radius. Right: The cumulative distribution of the virial parameter. Same colour scheme 
as used in the left-hand plot. 
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Figure 9. Fraction of MM-only sources as a function of the mass 
of the object. The error bars are identified by the grey boxes. 



dicating that the more massive clumps are slightly more 
gravitationally bounded. This translates as a slightly larger 
median virial parameter of 1 for the MM-only cores, com- 
pared with the star forming sources (median a — 0.9). 



6.3 MM-only sources in the context of massive 
star formation 

High-mass prestellar cores, or high-mass starless cores, the 
precursors of high-mass protostars in an analogous scaled-up 
version of low-mass starless cores, are still a missing piece of 
the high-mass star formation puzzle. Whether they may ex- 
ist will serve to constrain competing massive star formation 
theories which propose coalescence and accret ion as possible 
form ation mechanisms. Infrared dark clouds (jjackson et al.l 
2008) have been proposed as possible candidates to host such 



starless objects. However, high-angular resolution observa- 
tions of a small sample of candidate massive starless cores 
have syste matically revealed embedded, accreting protostars 
inside (e.g. iPillai et alj|2006h . 

We have identified a new class of source, the MM-only 
core, which are excellent candidates for early stage proto- 
stars or massive young stellar objects. SED analysis here 
reveals the MM-only sources to have a similar mass range, 
with a larger number of less massive sources, compared to 
the star forming sources, and very similar physical proper- 
ties as these sources. The similar distribution of the virial 
state of equilibrium suggests that the MM-only sources are 
just as likely to form stars. These results suggest that the 
MM-only sources may in fact be scaled down objects com- 
pared with the star forming sources. The smaller mass, radii 
and linewidths of the MM-only cores further suggests that 
they will possibly form intermediate mass stars rather than 
high-mass stars. 

Alternatively, the MM-only cores may not produce 
stars, and could remain starless transien t clumps (e.g. 
iBeuther fc Hennind[2009l ; iBuckle et al.ll2009T l. Further spec- 
tral line observations designed to ascertain the chemistry of 
the MM-only core are currently underway. 



6.4 Is the MM-only population comprised of 
sub-samples? 

It is also possible that the MM-only clumps could be com- 
prised of multiple sub-samples or populations, with one in- 
terpretation that each of these sub-populations likely follow 
different evolutionary paths. In this instance, they could be 
comprised of both cores with high and low mass stars, cores 
which are in the very earliest stages of evolution, and also 
possibly some cores that are not destined to begat high mass 
stars. 

Attempts to break class MM sources into multiple pop- 
ulations with respect to the mass, radius, linewidth, lumi- 
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nosity and temperature reveal the MM-only population to 
have a continuous distribution with respect to each of these 
parameters (see for instance Fig. [6] [7] and [8]). Consequently, 
there is currently no physical indication that the MM-only 
sources are comprised of distinct sub-populations. 

Despite this, from our results we can infer that the MM- 
only sources will have different star formation histories de- 
pending on their initial mass. Figure [9] presents a plot of the 
fraction of MM-only sources as a function of the mass of the 
clump. All of the sources in our sample less than lOOM©, 
are MM-only sources, and the fraction of objects that do 
not show any signature of star formation decreases as the 
mass of the clump increases. 

This indicates that the most massive clumps have a 
higher efficiency to form massive stars. Additionally it may 
suggest that there is a dichotomy in mass around 100- 
2OOM for the MM-only sample. The MM-only sources be- 
low this threshold will likely not form massive stars. 

Assuming an initial mass function (IMF) as per 
IChabrierl (|2003D . and a star formation efficiency of 50 
per cent, then we would need a minimum clump mass of 
~ 5OOM0 in order to support the formation of one star in 
excess of lOM© . That is, the maximum mass of the original 
clump dictates the size of stars that can form. Thus, clumps 
with masses lower than lOOM© will not form massive stars, 
whilst more massive MM-only sources on the contrary have 
enough material to form stars in excess of IOMq. 

Assuming that the evolution of protostellar objects is 
faster for more massive objects could explain why we detect 
a smaller fraction of MM-only sources i.e. before signatures 
of star formation, for increasingly massive objects. In this 
instance, the MM-only cores do represent targets of choice 
to assess the earliest stages of star formation. 



7 CONCLUSIONS 

We have undertaken an ammonia molecular line study, in 
the lowest two inversion transitions, of young massive star 
formation regions in the southern hemisphere as part of an 
effort to characterise the earliest stages of high-mass star 
formation. The sample targeted was derived from the mil - 
limetre continuum emission studies o f iHnT et al. I (|2005ll2006h 
and included sources with and without signatures of high 
mass star formation. In total, 244 sources were observed in 
both ammonia transitions, using the Parkes radio telescope. 
Of these, 138 had detections (>3-<r) in NH3(1,1), includ- 
ing two sources with two velocity components, and 102 in 
NH3(2,2). The MM-only sources are not more or less likely 
to be detected in either transition than the sources with a 
methanol maser and/or radio continuum association. 

The spectral line properties of the sources have been 
determined from gaussian fits to the lines: linewidth, flux 
density, opacity. In addition, physical parameters, such as 
the rotational and kinetic temperatures as well as the col- 
umn density, were derived from the spectra. From the ki- 
netic temp erature, and rev isiting previous SED modelling 
techniques l|Hill et al.ll2009h . we have determined the mass 
and luminosity of the sources in the sample. We have used 
the Bayesian inference method of SED modelling, which pro- 
vides robust estimates of the parameters as well as good es- 



timates of the uncertainties associated with each parameter, 
further allowing robust statistical conclusions. 

Combining continuum and line observations has proven 
to be quite powerful. In our case, ammonia observations and 
SED modelling are very complementary in terms of parame- 
ter determinations. It is clear, that the kinetic temperature, 
as derived from ammonia, used in combination with SED 
modelling has constrained the range of validity for the tem- 
perature, which in turn has constrained the range of validity 
for the mass and luminosity for each of the sources in our 
sample. 

The MM-only sources, those with no overt signs of 
massive star formation, have smaller NH3(1,1) and (2,2) 
linewidths on average compared with sources associated 
with a methanol maser and/or radio continuum source. They 
are also less massive and smaller on average but reach the 
same upper values for both parameters. Because the MM- 
only sources have similar brightness/flux, column densities 
and similar temperatures as star forming sources, they have 
the potential to form stars. The least massive MM-only 
sources, cannot form high mass stars and will likely proceed 
to form intermediate mass stars, whilst the more massive 
clumps, are more interestingly strong candidates for early 
stage massive protostars. The different hypotheses will re- 
sult in different internal structures. 

Higher spatial resolution observations with ALMA will 
allow resolution of the internal structures of the MM-only 
cores, which may then allow us to distinguish between star- 
less cores that will begat massive stars and those that are 
transient. 
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APPENDIX A: NH 3 (1,1) AND (2,2) SPECTRA 

The NHs(l,l) and (2,2) inversion spectra detected by the 
K-band receiver on the Parkes Telescope are presented here 
for all sources, including those in the sample figure (Fig. [IJ. 
Sources are presented in right ascension order, consistent 
with Tabl e [T]an d the complementary continuum images of 
iHill et all (|2005l ). The x-axis is in units of velocity (kms 1 ) 
whilst the y-axis is the flux density (Jy/beam) or strength of 
the detection. Both the Vlsr and flux density of the sources 
are presented in Table [T] The fits as returned from CLASS 
are overlaid on the spectra in green. 

This paper has been typeset from a Tp^X/ DTgX file prepared 
by the author. 



